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Abstract 


It  is  shown  that  the  small-angle  scattering  of  X-rays  or  neutrons  by  dislo- 
cations within  a deformed  metal,  which  are  partially  ordered  into  wall-like 
structures,  is  characterized  by  several  factors.  Principally  these  are  associated 
with:  1)  a single  dislocation  or  dipole,  2)  the  dislocation  configuration  in  the 
plane  of  the  wall,  and  3)  the  distribution  of  dislocations  across  the  wall  thick- 
ness. With  the  assumption  of  isotropic  elasticity,  small-angle  scattering  will 
be  sensitive  only  to  the  edge  components  of  the  dislocations.  The  scattered 
intensity  is  dominated  by  scattering  from  dislocations  which  he  perpendicular 
to  the  scattering  vector,  q,  and  reaches  a maximum  when  q is  normal  to  the 
slip  plane  of  these  dislocations.  Above  a particular  |q|,  the  scattered  intensity 
is  sensitive  only  to  the  total  edge  dislocation  content  of  the  visible  disloca- 
tions (i.e.  scattering  is  incoherent),  while  below  this  value,  the  scattering  is 
dominated  by  how  the  dislocations  are  distributed  in  walls.  For  walls  normal 
to  their  slip  planes,  the  configuration  factor  will  reflect  the  dislocation  distri- 
bution in  the  plane  of  the  wall,  while  for  walls  parallel  to  their  slip  planes,  the 
distribution  in  the  thickness  direction  will  be  visible.  Therefore,  even  though 
a deformed  material  is  composed  of  complicated  dislocation  structures,  only 
those  segments  conforming  to  these  rather  strict  prescriptions  will  be  singled 
out  for  scattering,  and  by  adjusting  the  beam/slip  system  geometry,  many 
parameters  of  the  microstructure  can  be  determined  experimentally. 
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I.  INTRODUCTION 


This  paper  is  motivated  by  experiments1  being  carried  out  on  the  NIST  materials  science 
beam  line  at  the  National  Synchrotron  Light  Source  at  Brookhaven  National  Laboratory, 
exploring  ultra-small-angle  X-ray  scattering  (USAXS)  by  deformed  single-crystal  metals. 
Here,  we  consider  the  theory  of  small-angle  scattering  from  individual  defects  (dislocations 
and  dislocation  dipoles)  and  present  a framework  for  understanding  the  role  of  wall  config- 
urations composed  of  these  defects. 

The  mechanical  response  of  a metal  to  plastic  deformation  has  traditionally  been  divided 
into  several  stages.  Stage  I is  characterized  by  a nearly  constant  flow  stress  as  a function  of 
applied  strain  and  is  only  observed  in  single  crystals  at  temperatures  much  lower  than  their 
melting  points.  In  stage  II,  the  stress/strain  behavior  becomes  linear  with  a significant  slope, 
resulting  in  work  hardening.  This  stage  II  slope  is  temperature  independent.  In  stage  III, 
temperature-dependent  recovery  processes  occur  resulting  in  a decreasing,  near-parabolic 
work  hardening.  The  existence  of  later  stages  of  work  hardening  have  also  been  reported, 
but  stages  II  and  III  are  the  most  important  industrially  and  discussion  is  restricted  to  these 
stages. 

According  to  the  current,  imperfect,  understanding  of  deformation  in  fee  metals,  the 
elastic  self  energy  of  the  system  drives  the  dislocations  to  organize  themselves  into  partially- 
ordered  wall  structures.  In  stage  II2-6,  where  our  interest  begins,  dislocation  walls  often 
develop  parallel  to  the  active  slip  planes,  perhaps  forming  partially-ordered  Taylor  lattices2. 
These  wall  structures  are  sometimes  referred  to  as  “carpets” . In  Stage  III,  the  walls  form  the 
boundaries  of  space-filling  cells2-6.  There  have  also  been  proposals,7,8  partially  supported 
by  observations,  that  the  walls  may  include  dislocation  dipoles.  Thus,  we  will  focus  on 
scattering  by  dislocations  distributed  in  walls  of  various  types. 

The  theory  of  small-angle  scattering  by  dislocations  has  been  the  subject  of  a number 
of  earlier  papers9-12,  stimulated  by  early  small-angle  scattering  experiments  with  X-rays 
(SAXS)  and  neutrons  (SANS).  However,  the  importance  of  the  wall  and  cellular  structures 
described  above  was  not  fully  appreciated  at  the  time  of  this  work,  so  the  theory  was 
necessarily  incomplete.  Seeger  and  Kroner9  carried  out  an  elegant  analysis  for  the  dilatation 
field  of  the  edge  dislocation  for  the  general  case  in  which  the  dislocation  is  not  straight. 
However,  since  we  will  be  working  with  mostly  2D  wall  structures,  we  find  the  2D  theory 
of  Atkinson  and  Hirsch12  more  in  tune  with  our  point  of  view,  and  we  will  refer  often  to 
it.  Certainly  dislocations  in  deformed  metals  are  not  straight.  Nevertheless,  2D  theory  is 
appropriate  for  this  problem  because  the  dislocations  roughly  approximate  straight  lines  over 
the  length  scales  probed  by  the  small-angle  scattering  experiments.  In  terms  of  the  theory 
presented  in  this  paper,  it  is  necessary  for  the  dislocations  to  be  approximately  straight  over 
length  scales  much  greater  than  the  inter-dislocation  distance. 

In  contrast  to  Bragg  peak  broadening13,  small-angle  scattering  is  caused  exclusively  by 
the  local  volume  change  induced  by  the  dislocations.  Thus,  in  the  linear  elastic  approx- 
imation, only  the  edge  components  of  dislocations  will  be  visible.  This  is  because  screw 
dislocations  have  no  net  volume  change.  We  focus  on  the  elastic  part  of  the  field  because 
the  additional  volume  changes  due  to  nonlinearity  in  the  strain  fields  of  both  edges  and 
screws  can  be  neglected  in  practice12  for  SAXS. 

In  the  kinematic  or  single-scattering  approximation,  the  scattered  intensity  is  composed 
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of  a product  of  several  factors:  a factor  for  the  individual  atoms,  one  for  a single  defect 
(either  single  dislocation  or  dislocation  dipole),  a factor  associated  with  a configuration  of 
identical  defects,  and,  finally,  a form  factor  for  the  entire  sample.  In  the  following  sections, 
we  address  each  of  these  in  turn,  except  for  the  atomic  factor,  which  in  the  case  of  a pure 
metal  can  not  provide  an  additional  scattering  contrast  mechanism,  and  the  sample  size 
factor,  which  appears  in  a scattering  vector  range  totally  unreachable  in  practice.  The 
analysis  sections  are  followed  by  a section  on  computer  simulations  of  relevant  dislocation 
structures  to  validate  the  analysis. 


II.  SINGLE  DEFECT  STRUCTURE  FACTORS 

The  scattering  at  small  angles  can  be  written  in  terms  of  the  total  scattering  amplitude 
function,  A,  as14 


/=  4(qM*(q), 

4(q)=  aaJ  (l  + <j)M (x))  exp(iq  ■ x)  dx,  (1) 

where  aa  is  the  scattering  amplitude  for  a unit  volume  of  perfect  crystal,  and  the  integral 
is  the  amplitude  structure  factor  for  the  scattering  volume.  The  dislocations  within  this 
volume  contribute  a total  fractional  density  change,  (f>tot,  at  the  point  x.  q is  the  scattering 
vector.  The  first  term  in  the  integrand  is  scattering  from  the  entire  sample  which  only 
appears  at  extremely  low  |q|.  It  does  not  contribute  to  the  visible  small-angle  scattering 
and  can  be  ignored.  Thus,  the  scattering  amplitude  becomes 

A(q)  = aa  J <j)tot{x)  exp(iq  • x)  dx.  (2) 

If  the  defects  are  distributed  over  the  volume  of  the  material  at  the  positions,  R*,  then 
the  fractional  density  change,  0tot(x)  is  the  linear  sum  of  the  contributions  from  each  defect, 
and  if  they  are  all  equivalent,  then  a decomposition  can  be  made, 

1 = -^  = \ad\2  c,  (3) 

laa| 


where 


ad=  J </>(x)  expjzq  • x}  dx 

c=  exP{*3  * (Ri  - Rj)}.  (4) 

Here,  ad  refers  to  the  structure  factor  of  a single  defect,  whether  a single  dislocation  or  a 
dislocation  pair,  etc.,  and  c is  a discrete  finite  sum  that  refers  to  the  distribution  of  defects. 
We  will  call  c the  “configuration  factor”  since  it  is  technically  not  a structure  factor.  <f>(x) 
is  the  fractional  density  change  (dilatation)  for  a single  defect,  with  the  defect  positioned  at 
the  origin.  Of  course,  in  a finite  sample,  some  of  the  defects  will  be  located  near  an  external 
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surface,  and  the  boundary  conditions  for  all  the  defects  are  thus  not  equivalent,  but  this 
effect  can  be  safely  ignored. 

This  decomposition  into  single-defect  and  configuration  factors  is  only  valid  if  all  the 
defects  are  the  same.  In  the  general  case,  dislocations  in  a given  configuration  will  possess 
different  Burgers  vectors  and  different  slip  planes,  so  they  are  not  equivalent.  In  a later 
section,  we  will  show  how  to  write  these  factors  so  that  they  can  be  decomposed  in  the 
desired  way;  for  now  we  work  out  the  form  for  aa,  as  defined  above,  to  see  how  it  depends 
on  Burgers  vector  and  slip  plane,  saving  the  generalization  for  later. 

The  first  case  to  be  treated  is  the  single  dislocation,  aligned  along  the  Z-axis.  For  most 
cases  of  small-angle  scattering,  it  can  be  assumed  that  the  scatterer  and  matrix  have  two 
distinct  densities  with  a sharp  boundary  between  them.  For  dislocations,  we  must  consider 
isotropic  elasticity,  where  only  the  edge  component  of  the  dislocation  contributes  a density 
change15,  given  by 


4>=  beK 


sin# 


k=  — 


u 

1 1 — 2u 
2n  1 — v ’ 


(5) 


where  the  edge  component  of  the  Burgers  vector,  6e,  is  along  the  X-axis,  and  the  polar 
coordinates  6 and  u specify  the  position  of  the  field  point  on  the  X-Y  plane  with  the  angle 
measured  with  respect  to  the  X-axis,  as  shown  in  Fig.  1.  u is  Poisson’s  ratio. 

As  noted  earlier,  the  dislocations  show  enough  two  dimensionality  to  make  2D  modeling 
a reasonable  approximation.  But  it  is  necessary  to  limit  the  distance  in  the  “Z-direction” . 
The  reason  is  not  that  a dislocation  ever  ends  in  the  crystal,  but  rather  that  a dislocation 
length  is  part  of  a large  loop  which  closes  on  itself  over  some  distance.  More  generally,  the 
dislocation  is  tied  into  a larger  structure  over  a distance  which  we  can  model  as  an  effective 
length  for  the  2D  dislocation.  Thus,  we  integrate  along  the  dislocation  line  (Z-direction)  for 
a finite  length,  — H to  H , 

r Tl  r2n  gjjj  Q rH 

o-singie=beK  / exp{iqpu  cos(9  — £)}tt  du  dd  / exp {iq2z}dz 

Jn0  Jo  xl  J—h 

A u . c \ Jo {Qp'R')  JoiQp'R'o) 

= 47r beK  sin  £ 

L Qp  Qp 

« 4 

Qp  Qz 

Here,  aSingie  is  the  structure  factor  for  a single  dislocation  and  J0  is  the  zero-order  Bessel 
function.  The  polar  angle  £ and  the  radial  distance  qp  give  the  projection  of  the  scattering 
vector  onto  the  X-Y  plane.  As  before,  the  angle  is  measured  with  respect  to  the  edge 
component  of  the  Burgers  vector  (X-axis).  In  these  equations,  the  last  step  is  valid  only 
when  qp  1Z  >•  1 and  qp TZq  <C  1,  where  1Z  is  the  upper  cut-off  for  the  dislocation  strain  field, 
and  IZo  is  the  lower  cut-off  at  the  effective  dislocation  core  radius.  For  this  derivation  to 
be  correct,  it  is  necessary  that  H TZ.  The  exact  expression  indicates  that  there  is  no 
scattering  by  the  dislocations  at  qp  = 0,  as  noted  by  Dexter16;  this  must  be  correct,  because 
the  dislocation  induces  zero  net  volume  change  in  the  material.  When  the  scattering  vector  is 


sin  (qzH) 

Qz 

(6) 
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very  small,  it  samples  over  a region  of  zero  volume  change,  and  is  not  scattered.  These  results 
are  in  the  form  already  given  by  Atkinson  and  Hirsch12.  We  note  here  the  very  important 
point  that  the  length  factor,  sin (qzH/)/qz,  is  a function  with  a very  strong  and  sharp  peak 
at  qz  = 0 for  dislocation  lengths  typical  of  deformed  metals.  Thus,  the  effective  dislocation 
length  can  be  directly  measured  in  experiments  that  probe  the  angular  dependence  of  the 
small-angle  scattering. 

The  second  case  to  be  treated  is  a pair  of  dislocations  displaced  from  one  another  in 
the  Y-direction  by  the  distance,  2d.  See  Fig.  2.  The  strain  field  of  the  two  dislocations  is 
expanded  in  a Taylor  expansion  up  to  first  order.  There  are  two  regions  of  expansion,  u > d 
and  u < d,  but  the  two  expansions  have  the  same  functional  form,  and  in  a straight  forward 
extension  of  the  above  analysis,  the  total  structure  factor  can  be  written 


® dipole  ~ 


8nbeK sin2 £ j2  sin (qzH) 
qp/2  V 7 id  qz 


(?) 


which  is  valid  in  the  same  range  of  g-space  as  before,  and  where  the  oscillatory  part  of  a 
Bessel  function  is  again  ignored. 

A comparison  of  Eqns.  (6)  and  (7)  shows  two  very  important  differences  between  scat- 
tering by  single  dislocations  and  scattering  by  dipoles.  First,  the  g-dependences  of  the 
scattering  amplitudes  differ  by  q~ !/2,  which  results  in  a q~l  difference  in  the  scattered  in- 
tensities. In  the  reciprocal  space  range  where  qp  « 2n/d  and  ignoring  angular  effects,  the 
dipole  and  the  single  dislocation  have  similar  scattering  intensities.  In  USAXS  experiments, 
however,  qp  <C  2tt /d,  thus  resulting  in  a much  higher  scattered  intensity  for  the  dipole.  This 
is  easily  understood  since  a dipole  can  be  regarded  as  a planar  region  of  vacancies  which  has 
a net  volume  change  at  small  gp,  whereas  the  dislocation  dilatation  field  averages  to  zero. 
The  requirement  that  qp  <C  27t /d  just  means  that  the  X-ray  beam  must  sample  a region 
larger  than  the  defect.  Second,  the  angular  distributions  of  the  scattering  are  profoundly 
different.  The  scattered  intensity  from  dipoles  goes  as  sin4£  as  compared  with  sin2£  for 
single  dislocations.  Thus,  the  high  relative  scattering  by  dipoles  over  single  dislocations  will 
only  be  visible  over  a narrow  range  of  orientational  angles. 

In  addition  to  the  elastic  dilatation,  which  is  a feature  of  only  the  edge  component  of  the 
dislocation  Burgers  vector,  dislocations  possess  density  changes  in  the  core,  and  the  nonlin- 
ear part  of  the  elastic  field  (even  the  shear  component)  also  contributes  a density  change. 
However,  as  has  been  shown  elsewhere,  these  contributions  are  negligible  compared  to  the 
elastic  expressions  in  the  region  of  qp  space  accessible  to  ultra-small-angle  scattering.12  The 
same  is  true  of  vacancies  and  small  collections  of  vacancies  that  can  form  during  deformation. 


III.  CONFIGURATION  STRUCTURE  FACTORS 
A.  General 

The  main  difference  between  this  work  and  that  of  previous  authors  is  that  we  use  a dif- 
ferent factor  for  the  dislocation  configuration.  Atkinson  and  Hirsch12  treated  the  dislocations 
as  a cloud  of  independent  scatterers.  Seeger10  treated  dislocation  pile-ups,  but  these  config- 
urations are  rarely  observed  in  deformed  metals.  Seeger  and  Brand11  also  considered  small 
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dislocation  loops  scattering  independently,  which  is  appropriate  for  an  irradiated  sample.  In 
fact,  it  is  now  well  understood  that  dislocations  in  deformed  metals  form  partially-ordered 
structures,  mostly  walls,  and  that  is  the  case  considered  here. 

Dislocations  in  a wall  can  be  divided  into  distinct  populations.  First,  for  a given  slip 
system,  there  exist  dislocations  with  both  positive  and  negative  Burgers  vectors.  Second, 
dislocations  from  more  than  one  slip  system  may  be  present  in  a wall,  depending  on  the 
deformation  history  and  the  strain  geometry.  Third,  the  dislocations  may  lie  in  different 
directions  in  the  wall.  Finally,  some  of  the  dislocations  may  be  better  represented  as  dislo- 
cation dipoles,  rather  than  single  dislocations.  The  analysis  for  this  completely  general  case 
is  quite  complicated,  and  simplifying  assumptions  must  be  made  to  build  a useful  physical 
picture.  In  particular,  to  obtain  a usable  solution,  it  is  necessary  to  achieve  a decomposition 
into  a number  of  factors,  and  the  first  goal  in  this  section  is  to  show  under  what  conditions 
this  can  be  achieved  for  the  single  defect  and  configuration  factors. 

Factorization  cannot  be  achieved  if  both  dipoles  and  dislocations  coexist  in  a wall,  be- 
cause the  single  defect  functional  dependence  on  qp  is  different  for  the  two  cases.  But  if  only 
one  class  of  defect,  either  single  dislocations  or  dipoles,  is  dominant,  then  the  functional 
form  for  this  scattering  can  be  separated  from  all  other  configurational  considerations.  As 
shall  be  seen  later,  however,  there  is  a range  in  reciprocal  space  where  the  contributions  of 
single  dislocations  and  dipoles  can  be  separated  and  direct  experimental  measurements  of 
the  relative  numbers  of  these  defects  may  be  possible. 

The  problem  can  be  simplified  further  by  focusing  on  single-crystal  fee  metals,  with 
{111} (110)  slip  systems.  In  the  early  stages  of  deformation,  i.e.  Stage  II,  walls  often 
develop  parallel  to  the  {111}  slip  planes — the  “carpets”  referred  to  in  the  Introduction. 
The  dislocations  are  generally  parallel  to  one  another,  and  may  take  the  form  of  partially- 
ordered  Taylor  lattices2.  In  later  stages  of  deformation,  walls  are  formed  roughly  normal 
to  the  primary  slip  planes,  presumably  the  end  result  of  cross  slip  and  climb  from  the 
initial  slip  planes,  or  the  interaction  between  two  slip  planes  operating  simultaneously.  For 
example,  two  separate  slip  planes  will  intersect  in  (110)  type  directions.  Dislocation  walls 
created  by  such  an  interaction  contain  Burgers  vectors  from  two  different  slip  planes,  but 
the  dislocations  will  all  be  parallel  to  the  intersection  line  of  the  two  slip  planes.  Thus,  for 
either  the  carpets  or  the  stage  III  walls,  it  is  reasonable  to  assume  that  the  dislocations  in  a 
given  wall  all  lie  essentially  parallel  to  one  another,  even  though  they  will  have  both  Burgers 
vector  sign,  and  may  even  have  more  than  one  Burgers  vector  direction  on  more  than  one 
slip  plane.  If  the  defects  all  lie  parallel  to  one  another,  it  is  possible  to  factor  out  the  line 
direction  part  of  the  scattering  function,  and  a further  factorization  has  been  achieved. 

We  start  by  rewriting  the  structure  factor  for  an  individual  defect,  i,  as  a product, 

Of  = PiAiMO-  (8) 

In  this  equation,  since  the  defects  can  belong  to  more  than  one  slip  system,  they  may  have 
a variety  of  Burgers  vector  edge  components,  Pi  = ±(6e)i,  and  angular  factors,  /*(£).  As 
shown  in  Eqns.  (6,  7),  / and  A have  different  functional  forms  for  single  dislocations  and 
for  dipoles.  However,  A does  not  depend  on  the  sign  of  the  Burgers  vector  or  the  particular 
slip  plane.  With  the  above  notation,  we  include  the  Burgers  vector  and  angular  factors,  Pi 
and  /;(£),  in  the  configuration  factors,  and  write  W for  the  new  configuration  factor, 


6 


(9) 


W = Y PiPjfi(€)fj(t)  exp{«q  • (Ri  - Rj)}, 
id 


where,  from  Eqns.  (6,7), 


1=  A2W 

47 tk  f sin (qzH) 

- single ~ I 

9P  V 92 

87t k / 2 / sin(< 

9p/2  V 7T0?  V 9 

Although  this  separation  of  A and  W is  only  rigorously  valid  for  materials  that  satisfy 
the  earlier  assumption  of  isotropic  elasticity,  the  general  results  of  this  paper  should  be 
substantially  the  same  for  anisotropic  crystals. 

These  equations  are  valid  only  if  a particular  defect  type  (dipole  or  single  dislocation)  is 
dominant  in  a given  wall.  If  both  dipoles  and  single  dislocations  exist  in  significant  numbers 
in  the  same  wall,  then  the  scattered  intensity  cannot  be  factored  in  this  way,  and  interference 
between  the  two  types  of  defects  will  be  present.  A second  proviso  is  that  Eqn.  (9)  applies 
to  a single  wall.  Scattered  intensity  from  more  than  one  wall  is  assumed  to  be  incoherent, 
and  is  simply  a sum  over  the  different  walls.  For  q < 0.0002A-1,  at  large  strains  in  Cu 
or  A1  solid  solution  alloys,  complete  cells  will  scatter;  but  this  type  of  scattering  is  not  yet 
experimentally  accessible,  and  is  not  the  subject  of  this  paper. 

Equations  (9)  and  (10)  describe  the  scattering  from  a specific  configuration  of  dislocations 
or  dipoles  in  a particular  wall.  For  the  averaged  scattering  over  a collection  of  similar  walls, 
the  averaged  intensity  becomes 


(I)  = A2  (W) 


(11) 


Now  it  is  clear  why  the  defect  structure  factor,  a^,  introduced  in  Eqn.  (3),  is  not  the 
appropriate  expression  when  averages  over  probability  distributions  must  be  taken;  only  A 
is  invariant  for  all  defects  of  a given  type  (i.e.  single  dislocation  or  dipole)  in  a distribution 
of  defects. 

The  analysis  of  (W)  forms  the  content  of  the  remainder  of  this  subsection.  First,  a joint 
pair  distribution  function  is  defined  for  the  double  sum  over  pairs  in  Eqn.  (9).  Let  p(R*,  Rj) 
be  defined  as  the  joint  probability  for  finding  one  dislocation,  labeled  i at  Rj  and  a second 
dislocation  labeled  j at  Rj.  Then,  following  (9),  (W)  is  written 


<w>  = Y 0i0jfi(£)fj(€)  j J p(R*,  Rj)  exp{iq  • (Rj  - Rj)}  dR,  dRj 

id 

J J p(Ri,  Rj)  dKi  dRj  = 1,  (12) 


where  the  second  equation  gives  the  normalization  used  for  the  joint  pair  distribution  func- 
tion. In  the  double  sum  over  the  dislocation  pairs,  the  diagonal  terms  for  i = j refer  to  the 
same  dislocation,  and  represent  the  self  scattering  of  the  separate  dislocations.  For  these 
self  scattering  terms,  of  course,  p = <5(Ri,  Rj).  For  all  pairs,  however,  where  the  dislocations 
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making  up  the  pair  are  different  dislocations,  the  joint  pair  distribution  function  will  be 
assumed  to  be  a smooth  function.  Thus,  we  can  write 


+ Y/0i0jfM)fj(O  [ /p(Ri, R,)exp{iq- (R,  -R,)}dRjdR,.  (13) 


Now  remember  that  the  dislocations  consist  of  several  populations  with  different  Burgers 
vectors,  slip  planes,  etc.,  and  that  within  each  type,  e.g.  a given  Burgers  vector,  etc.,  the 
dislocations  are  equivalent  and  have  equal  p(R*,Rj).  Thus  the  sums  over  dislocations  of  a 
given  type  can  be  replaced  by  Nm,  etc.,  where  Nm  refers  to  the  number  of  dislocations  of 
type  m.  Accordingly,  the  component  parts  of  (W)  = (Wmm)  + (W^),  can  be  written, 

<wmm)  = y NmfmfKo 

771 

+ J2Nrn(Nm-l)Pmfm(€)  / [ R-U  exp{«q  • (R^  - R^)}  dRm  dR'm, 

771  ^ ^ 

<Wm„>=  y NmNA/3nfm(()fn( f) 

m^n 

x J J V{R m,  R4)  exp{?q  • (R™  - Rj,)}  (C  dR'n,  (14) 

where  the  first  equation  refers  to  dislocations  belonging  to  the  same  type,  and  the  second 
to  dislocations  of  different  types.  The  factor,  Nm(Nm  — 1)  in  this  equation  arises  because 
there  are  Nm(Nm  — 1)  terms  in  the  off-diagonal  double  sum  for  each  dislocation  type  in 
the  second  term  of  (13).  The  new  joint  pair  distribution  function,  ^(R^R^)  is  now  the 
joint  probability  of  finding  one  dislocation  of  type  m at  position  R™  and  a second  (different) 
dislocation  of  type  n at  R^.  (Of  course,  the  dislocations  may  be  of  the  same  type,  m = n 
or  of  different  types,  m ± n.)  It  is  very  important  to  note  that  this  joint  pair  distribution 
function  is  not  the  standard  pair  correlation  function.  These  equations  can  be  combined 
into  the  final  equations, 


<W)  = £ *«/£/£(€)  [l  - / f P(Rm,  RJJ  exp{iq  ■ (R„  - 

771  l J J J 

+ NmNnpm0nfm (€)fn (f ) f [ V(K m,  RJJ  expjzq  • (R^  - Rjj}  dR’n, 

m,n  ** 

J J ^(Rm,R^)  dKmdR'n  = 1, 


(15) 


where  the  double  sum  is  taken  over  all  (m,  n) . The  last  equation  sets  the  normalization  for 
the  joint  pair  distribution  functions. 

Eqn.  (15)  has  a very  important  form.  The  first  term,  the  one  with  unity  in  the  brackets, 
is  the  self  scattering  term,  which  is  independent  of  the  distribution  function  because  a given 
dislocation  is  fully  correlated  with  itself.  It  is  also  independent  of  the  scattering  vector,  q, 
and  is  proportional  to  the  number  of  scatterers,  Nm.  This  first  term  will  be  seen  to  have 
special  importance  in  the  following  for  the  case  where  the  net  Burgers  vector  in  a wall  is 
near  zero. 
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Finally,  we  note  that  although  it  is  physically  impossible  for  two  different  dislocations 
to  occupy  the  same  position,  this  does  not  mean  that  we  must  require  T^R^R^)  = 0 
for  Rm  = R4.  Although  physically  incorrect,  we  will  show  in  a later  section  that  the 
resulting  error  is  insignificant.  This  allows  us  to  use  much  simpler  distribution  functions  in 
the  analysis. 


B.  Walls  of  Zero  Thickness 


Because  the  R dependence  of  the  exponential  in  Eqn.  (15)  is  only  a function  of  the 
difference  between  the  dislocation  positions,  it  is  natural  to  introduce  a change  of  variables, 


Rm  + Rn~~  ^ rmn 

FLm  ~ R'n=  rmn-  (16) 


Here,  it  is  assumed  that  the  wall  has  zero  thickness;  in  the  next  subsection,  the  results  are 
generalized  to  finite  thickness.  Thus,  the  position  vectors  in  Eqn.  (15)  become  scalars,  and 
the  product  in  the  exponential  will  be  chosen  such  that  qw  lies  in  the  wall,  normal  to  the 
dislocation  lines.  Because  the  integration  is  finite,  one  must  integrate  over  the  appropriate 
limits.  For  a wall  which  extends  from  R = 0 to  R = +L,  the  initial  and  transformed 
integration  spaces  are  shown  in  Fig.  3.  Note  that  the  Jacobian  for  the  transformation  is 
1,  and  dRdR!  = dr  dr' . It  is  reasonable  to  assume  that  the  probability  distribution  for  the 
dislocations  depends  only  on  their  relative  positions,  r,  and  not  on  their  average  positions, 
r' . With  this  assumption,  the  integral  over  r'  is 

IIV Rn)  COS {qw(Rm  Rn)}  dRm  dRn 

= J J ^mn 

= fv(r  mn  ) cos {qyjrmn}  drmn  J drmn.  (1 0 

The  two  integrals  in  the  last  line  are  to  be  integrated  over  the  rotated  rhombus  in  Fig.  3. 
Because  of  the  symmetry  about  the  r'mn  axis,  the  integral  only  needs  to  be  carried  out  over 
the  upper  triangle  of  the  rhombus,  for  rmn  > 0.  Thus, 

fv(  t mn)  cos {qwr  mn  } dr 

mn  / cR  mn 

rL  r L rmn/2 

= 2 / V{rmn)  cos{qwrmn } drmn  / drmn 

JO  Jrmn/2 

rL 

= 2 / (L  ^mn)^(^ran)  COS-(^ti,7'TOn|-  drmn.  (16) 

Jo 

The  normalization  of  V(rmn ) follows  from  that  for  the  total  distribution  function, 


1=||  P(Rm,R'n)dRmdK  = J jv(  r mn  •>  'C mn  ) R mn  dr mn 

= 2 [ V{rmn){L  - rmn)  drmn ■ (19) 

Jo 
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When  all  dislocations  are  completely  random,  it  follows  that  VTandvm  — 1 /T2. 

The  function,  V(rmn),  requires  special  attention.  From  its  definition,  V(Rm,R'n ) = 
'P{rmn,r,mn)  = 'P(rmn)i  it  remains  a joint  pair  distribution  function,  even  though  it  depends 
only  on  the  single  variable,  rmn,  the  pair  separation  distance.  Thus,  V^mn)  should  not  be 
interpreted  as  the  probability  of  finding  a dislocation  pair  with  separation  rmn. 

Using  Eqns.  (17)  and  (18)  in  Eqn.  (15),  the  configuration  factor,  (W),  becomes 

<W>  = Nm0mfm(£)  (1  - Wnm»  + £ NmNn0m/3nfm(i)fn(t)  <Wm„> 

m m,n 

fL 

(Wmn)  = 2 (. L - rmn)V{  T mn  ) cos(qwrmn)  drmn.  (20) 

Jo 

This  is  a very  general  equation,  depending  only  upon  the  joint  pair  distribution  functions  for 
the  pairs,  (mn),  within  the  populations.  Note  that,  in  general,  correlations  exist  between 
pairs  of  plus  and  minus  Burgers  vectors,  and  between  dislocation  pairs  on  different  slip 
planes.  Recall  also  that  it  is  assumed  that  only  the  distances  between  the  dislocations  of  a 
pair  are  correlated.  The  midpoint  of  the  pair  is  randomly  distributed  over  the  wall,  and  is 
uncorrelated  with  other  pair  midpoints. 

The  part  of  the  sum  over  populations  which  relates  to  Burgers  vector  sign  can  be  further 
manipulated  if  the  excess  sign  (proportional  to  the  net  Burgers  vector),  Bs,  and  the  total 
number  of  dislocations,  Ns  for  each  slip  system,  s,  are  introduced, 

B,=  N*  - iv; 

NS=N+  + N-.  (21) 


Then  the  sum  over  the  positive  and  negative  dislocations  for  a given  slip  plane  group  can 
be  carried  out  explicitly  to  give 


<W>  = £ 62(s)/2«)  (N,  - JV+  (w+,+)  - N7  <W--}) 


+ £ 


Ms)W«U(C)/.'(f) 


(N,N,  + BA)  ((w;/)  + (w-r)) 

+ (NA  + ns,bs)  ((wsy)  - (w-r))  - 2 (NA  - BA)  (w+-) 

(Ws+/}  = 2 / (L  - r,,')P++ (r„i)  cos(qwrss’)  dr,,’  ,ete., 


(22) 


where  the  configuration  factors  { t "j  refer  to  both  dislocations  positive  and  'j  refer 
to  unlike  sign  Burgers  vectors.  It  is  assumed  that  the  configuration  factors  satisfy  the 
relations 


(wtf)  = (w-/) . (23) 

Again,  these  equations  are  quite  general,  and  describe  different  distributions  of  like  and 
unlike  signs  within  a given  slip  plane  set,  and  of  course  also  allow  for  different  distribu- 
tions between  slip  planes.  The  equations  are  complicated,  but  that  is  a reflection  of  the 
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complicated  nature  of  the  dislocation  ordering  problem.  For  example,  when  one  sign  of 
Burgers  vector  (say  — ) is  in  the  minority,  then  ordering  will  be  dominated  by  the  majority 
population,  and  the  minority  will  tend  to  be  disordered  (relative  to  itself),  with  possibly 
strong  correlations  between  unlike  dislocations.  Thus,  (W++)  and  (W+-)  may  be  strongly 
correlated,  while  (W~~)  may  be  nearly  random  since  the  average  negative  dislocation  sepa- 
ration greatly  exceeds  the  relevant  correlation  distances.  On  the  other  hand,  for  B = 0,  all 
distribution  functions  are  likely  to  be  strongly  correlated.  Since  this  B = 0 case  is  of  great 
interest  for  deformed  metals,  we  return  to  it  at  the  end  of  §IIID. 

When  the  populations  are  fully  disordered,  then  a good  approximation  will  be 

<VW?,)  = (ws+/)  = (w-r)  = (w+7 ) , (24) 

and  one  can  write  the  very  simple  result 

<w>  = £ bl(s)Nsfs(£)  (l  - (W?,))  + £ be(s)b'(J)f,(S)fX()B,B,  ,)  . (25) 

s s,s’ 

Thus,  if  there  are  no  special  correlations  between  like  and  unlike  sign  dislocations,  the  walls 
are  most  transparent  when  the  net  number  of  positive  and  negative  dislocations  are  equal 
and  Bs  = Bs>  = 0,  with  only  the  term  linear  in  the  N contributing.  This  linear  term  can  be 
simplified  further  by  solving  explicitly  for  (VV°S).  For  completely  uncorrelated  dislocations 
within  a single  slip  plane, 

/ o \ = sin2 (qwL/ 2) 

V ss/  qlL2/4  ’ 

<W)  = blNf(S)  (l  - ' (26) 

This  equation  has  some  very  interesting  properties.  First,  as  qw  — » 0,  the  scattering  goes 
to  zero  quadratically.  Second,  as  the  scattering  vector  increases  such  that  qwL  1,  the 
scattering  goes  to  the  self  scattering  limit,  proportional  to  the  total  number  of  dislocations, 
N.  Experimentally,  it  will  usually  be  the  case  that  qwL  >>  1,  so  the  wall  function  for  B — 0 
is  a constant,  independent  of  qw.  The  fact  that  the  wall  function  goes  to  zero  at  qw  = 0 
in  this  case  is  simply  due  to  the  fact  that  we  are  then  sampling  over  a region  with  no  net 
volume  change.  The  reader  should  not  confuse  this  behavior  of  (W)  as  qw  — >•  0 with  the 
near  zero  behavior  of  the  single  dislocation  structure  factor  of  §11.  In  particular,  qw  in  Eqn. 
(26)  and  qp  in  Eqn.  (6)  are  generally  not  equal  and  the  widths  of  the  two  holes  at  the  origin 
are  also  different. 

During  the  early  stages  of  uniaxial  deformation,  the  walls  contain  nearly  equal  numbers 
of  positive  and  negative  dislocations.  But  one  does  not  expect  the  configuration  factors 
to  be  exactly  equal  for  the  (++)  and  (-1 — ) cases.  We  return  later  to  the  case  where  the 
dislocations  are  partially  ordered. 

Now  compare  this  result  for  B = 0 with  the  case  of  a totally  disordered  wall  containing 
only  positive  dislocations  from  a single  slip-system  ( B = N).  Equation  (25)  can  then  be 
integrated  immediately,  with  V(r)  = 1/L2  to  give, 

l2 a, /2/v\  (i  sin 2(9wL/2)\  ,2„2f2K,  sin2(9„i/2) 

W (?)  l1  " llLVi  ) + b‘N  1 (fl  4^/4  ' 
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The  second  term  only  dominates  if  qw  <C  N/L 2 . For  realistic  experimental  conditions,  this 
requirement  can  only  be  satisfied  if  qw  is  nearly  perpendicular  to  q.  The  second  term  is  the 
standard  optics  diffraction  function  for  a single  slit  the  size  of  the  wall. 

Thus,  when  the  dislocations  are  totally  random,  and  qw  \JN/L 2,  the  scattering  is  the 
same  for  both  B = 0 and  for  B = N.  As  qw  —>  0,  the  scattering  will  be  proportional  to  the 
square  of  the  uncompensated  Burgers  vector. 

In  using  these  equations,  one  should  be  guided  by  the  apparent  fact  that  in  a deforming 
metal,  any  given  region  will  tend  to  be  in  single  slip,  with  a preponderance  of  dislocations 
from  a single  slip  plane.8  The  conjugate  systems  are  probably  important  as  a means  to 
form  “hard  spots”  which  tie  the  configuration  down  to  the  underlying  lattice,  but  usually  in 
any  given  local  region,  these  dislocations  are  only  a small  minority  in  the  total  dislocation 
population.  Thus,  for  scattering  purposes,  it  may  be  appropriate  to  use  the  equations  as 
if  only  a single  slip  system  is  operating  for  a given  wall,  and  average  over  the  different  slip 
systems  only  when  the  intensities  of  individual  walls  are  averaged. 


C.  Walls  of  Non-zero  Thickness 


The  walls  of  dislocations  formed  during  deformation  are  probably  diffuse  and  have  a 
nonzero  thickness.  In  general,  the  mechanisms  responsible  for  wall  formation  are  different 
in  the  length  and  thickness  directions.  For  example,  if  a wall  is  perpendicular  to  the  slip 
plane  of  the  component  dislocations,  then  the  formation  mechanism  in  the  length  direction 
would  include  dislocation  climb  and  source  distributions  whereas  the  thickness  formation 
mechanism  would  include  elastic  interactions  and  glide.  Therefore,  it  is  reasonable  to  assume 
that  the  dislocation  distributions  in  the  thickness  direction  and  in  the  length  direction  are 
independent,  so  that  V(w,t;w' ,t!)  = V(w,w')V(t,tr),  where  t and  w are  the  coordinates  of 
a dislocation  in  the  thickness  and  length  directions,  respectively.  Then  (W)  can  be  written 


<w>  = Y.  N<nPmfm(()  “ 2 / / V(Wm,  ®m)  COs(q,Xmm)  dwm  dw’m 

xf(L-ymm)V(y  mm  ) cos  {q^ym  m ) dymn^j 

T2  y ] NmNnPrnftnfrn(£)fn{£)  j j B(wm,  Wn ) COS{qtXmn ) dwm  dwn 
fL 

X / (L  ymnfB  {Umn  ) cos(qwymn)  dymn , 

Jo 


(28) 


where 


IV m VJn — Xmn 

+ VJn  2xmn 
tm  ~ Umn 

tm  + t’n=  2y’mn-  (29) 

are  the  equations  for  transforming  the  coordinates  into  the  “rotated”  coordinates  suitable 
for  the  integrations,  just  as  before.  To  simplify  the  analysis,  we  split  the  total  wall  factor 
into  factors  for  the  thickness,  (W^n),  and  length,  (W^n),  components, 
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<w>  = £*„££«)  (i  - (>*L,)  mJ) 

m 


+ E NmNm?Mm(OUO  (WL>  <vc,> 

tJijTl 


(wL)  - f f 'P(wm,w'n)cos (qtxmn)dwmdw’n 
(VO  = 2 [ (L-ym„)V(y  mn  ) cos{qwymn)  dy 

r 

J 0 


(30) 


We  will  explore  two  simple  cases  for  (W^n).  In  the  first  case,  the  dislocations  are 
randomly  distributed  over  a fixed  width.  Then  V(wm,w'n)  — \/l\,  for  0 < wm,w'n  < lt, 
where  /f  is  the  wall  thickness.  Such  a wall  has  sharp  boundaries  at  w = 0,  lt.  Then,  when 
the  coordinate  transformation  is  invoked  for  wm,  w'n,  as  before,  and  the  integrations  are 
carried  out  over  xmn  and  x'mn, 


sin2(qtlt/2) 


(31) 


This  thickness  factor  is  the  familiar  single-slit  diffraction  function  that  was  found  in  Eqn. 
(26),  and  has  a power  law  fall  off  at  high  qt  of  q t~2.  If  the  wall  is  sufficiently  thin  such  that 
qtlt  1,  then  (VV4n)  = 1,  and  the  result  reverts  to  that  of  the  previous  section. 

A second  simple  case  for  a wall  of  finite  thickness  can  be  obtained  by  assuming  that  the 
distribution  function  in  the  width  direction  is  Gaussian.  Thus, 


V(wm)  = (32) 

ifV7r 

where  lt  is  the  width  of  the  distribution,  and  the  total  joint  pair  distribution  function  in  the 
thickness  direction  is  V(wm)V(w'n).  Taking  the  same  transformation  of  coordinates  wm,w'n 
into  xmn,x'mn  as  given  in  Eqn.  (29),  the  joint  pair  distribution  in  xmn,x'mn  space  is 


V(x 


mn  i %mn 


i)  = V{wm)V{w'n)  = —2  exp 

7T  lt 


x 

mn  1 ^ mn 


/2) 


(33) 


so  that  the  width  part  of  the  distribution  function  is  again  a product,  and  (W^n)  can  be 
written 


. . r co  r oo 

(y^mn) — / / Xmn)  COs(.qt%mn)  dxmn  dxn 

' ' J —oo  J —oo 


cos (qtXmn)  dx  mn 


= \/2exp  ■ 


(34) 


This  equation  is  similar  to  the  sharp  (thick)  wall  case,  except  that  (W^n)  is  a Gaussian 
instead  of  a single-slit  diffraction  function.  Both  distributions  have  widths  of  lt,  but  the 
scattering  by  the  sharp  wall  has  a 1 /q?  power  law  fall  off  at  large  qt,  while  the  scattering  by 
the  diffuse  wall  decreases  more  rapidly. 
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It  is  important  to  note  an  essential  constraint  in  the  analysis  of  (W4n}.  In  the  earlier 
analysis  of  (W®n),  we  made  the  assumption  that  the  distribution  function,  V(ymn,y'mn)i 
can  be  factored,  'P(ymn,y,mn)  = ^>{ymn)'P{y'mn).  This  means  that  there  is  no  correlation 
between  the  distribution  of  pair  distances  and  pair  centers.  But  when  a pair  is  near  the 
end  of  the  wall,  this  assumption  must  break  down,  because  a pair  must  sense  the  end  of  the 
wall,  and  must  not  overlap  it.  That  is,  the  joint  pair  distribution  function,  V(ymn)  must  also 
depend  on  the  position  of  the  midpoint,  y'mn.  This  error  is  not  important  if,  as  expected, 
the  length  of  the  wall  is  very  large  compared  to  the  dislocation  correlation  length.  In  this 
case,  only  a very  small  percentage  of  the  dislocations  fail  to  follow  the  factored  distribution 
function.  But  for  the  wall  width  analysis,  the  wall  width  may  not  be  very  much  larger 
than  the  corresponding  dislocation  correlation  length.  In  this  case,  it  is  not  valid  to  split 
V(xmn,  x'mn)  into  independent  factors.  Thus,  except  for  special  cases  like  those  explored 
above,  one  cannot  go  further  than  the  first  equation  of  the  width  analysis,  Eqn.  (28),  in  a 
rigorous  general  treatment. 

Another  approach  for  the  treatment19  of  diffuse  boundaries  in  small-angle  scattering  has 
proven  very  useful  for  the  interpretation  of  scattering  by  block  co-polymers  with  diffuse 
interfaces.  We  now  adapt  that  approach  to  the  case  of  scattering  by  a dislocation  wall. 

If  the  dislocation  pair  distribution  function  is  uncorrelated,  it  follows  that  V(wm,w')  = 
V{wrn)V{w'n ),  and 


(win)  = / / 'P{wm,  w'n)  cos (qtxmn)  dwm  dw'n 


= J [f'PiXmn  + w'n)V{w'n)  dw'n 
= T{V®V). 


COS dx  r 


(35) 


In  this  equation,  T is  the  (ID)  Fourier  transform,  and  0 denotes  the  Fourier  convolution 
of  two  functions.  Following  Ruland20  (see  also  Koberstein,  et  al.19),  a wall  with  a diffuse 
boundary  is  introduced  as  the  convolution  between  a probability  function,  'Po(wm),  rep- 
resenting a sharp  wall  with  width  lt , and  a smoothing  function,  S(wm).  The  smoothing 
function  goes  uniformly  to  zero  for  large  wm  and  has  a diffuse  width,  ltd ■ The  probability 
function  for  the  non-sharp  wall  thus  has  the  form 


V{wm)  = S(wm ) 0 V0{Wm). 


(36) 


Using  the  Fourier  convolution  theorem,  we  obtain 

(™L)  = F(-P°  ® 5®  7>o  ® S)  = (WZ)  [^(5)] 2 * [^(S)]2,  (37) 

where  (VV^n)  is  the  sharp  wall  factor.  The  last  expression  is  valid  for  large  qtlt,  and  [^r(5')| 
will  be  a function  of  qt  which  sharpens  the  q~2  law  for  large  qt. 

Given  a model  for  the  smoothing  function,  5,  the  observed  deviation  from  q~2  scattering 
can  be  used  to  obtain  an  estimate  for  ltd-  This  issue  has  been  addressed  extensively,  most 
notably  again  by  Ruland20  and  Koberstein,  et  al.19.  The  most  common  approach  is  to  use 
the  sigmoidal-gradient  model  in  which  the  smoothing  function  is  assumed  to  be  Gaussian. 
The  resulting  thickness  part  of  the  wall  factor  is  then  given  by 
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where  a is  the  standard  deviation  of  the  Gaussian  smoothing  function. 

In  summary,  we  have  shown  that  reasonable  assumptions  about  the  joint  pair  distribution 
functions  in  walls  with  non-zero  thickness  lead  to  a configuration  factor  that  contains  a 
“width  factor”  in  addition  to  the  “length  factor”  discussed  in  the  previous  subsection.  For 
walls  with  sharp  boundaries,  this  width  factor  has  a q t-2  dependence  in  the  large  qt  (Porod) 
regime.  If  the  wall  has  a diffuse  boundary,  there  is  an  additional  qt  dependence  which  can  be 
quantified  in  terms  of  the  smoothing  function,  S(w).  This  additional  qt  dependence  serves  to 
increase  the  absolute  magnitude  of  the  slope  of  the  scattering  curve  in  g-space.  The  Gaussian 
case,  which  was  worked  out  in  detail,  corresponds  to  a Gaussian  smoothing  function  with  a 
wall  of  zero  width.  All  of  these  findings  are  in  keeping  with  classical  small-angle  scattering. 
However,  to  reach  these  results,  it  was  necessary  to  postulate  zero  correlation  between  the 
dislocations  in  the  width  direction.  This  may  be  a reasonable  assumption  in  some  cases, 
and  where  it  breaks  down  drastically,  it  may  be  possible  to  revert  to  the  dipole  treatment 

of  §11. 


D.  Ordered  and  Partially  Ordered  Walls 


The  long  range  interaction  of  the  dislocations  not  only  causes  the  dislocations  to  form 
wall  structures,  but  also  produces  varying  degrees  of  dislocation  ordering  within  the  walls 
themselves.  Such  internally-ordered  walls  scatter  differently  from  disordered  walls.  In  this 
subsection,  some  generic  effects  of  such  ordering  within  the  walls  will  be  explored. 

Consider  first  the  case  of  a totally  ordered  tilt  wall,  with  N+l  dislocations  of  the  same 
sign,  separated  from  one  another  by  the  distance,  L/N, 


= blf\0' 


N 


Y exp{iqwnL/N} 

,71=0 

sin2(qw(N-hl)L/2N) 


' N 

Y exp {-iqwnL/N} 

.71=0 


b2N2f2(£) 


sin  2{qwL/2N) 
sin2  (qwL/ 2) 
qIL2/ 4 ’ 


{N  > qwL ), 


(39) 


which  is  the  scattering  function  for  a grating.  For  qw  -C  N/L,  this  result  becomes  indis- 
tinguishable from  the  case  of  a random  wall  with  B = N,  and  qw  <C  ^Jn/L2.  See  Eqn. 
(27). 

A completely  ordered  B = 0 wall  with  alternating  sign  dislocations  gives 
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J^(-l)nexp  {iqwnL/N} 
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(40) 


Examination  of  Eqns.  (39)  and  (40)  shows  that  the  effect  of  replacing  every  other  positive 
dislocation  with  a negative  dislocation  is  to  “shift”  the  diffraction  pattern  by  half  the  period, 
or  7tN/L.  The  reason  for  this  unusual  behavior  is  that  each  positive/negative  dislocation 
pair  forms  a double-sized  unit  cell  with  peaks  spaced  half  as  far  apart  in  reciprocal  space. 
The  multiplying  form  factor  of  the  dislocation  pair  has  zeros  at  the  original  lattice  sites, 
causing  the  apparent  shift.  The  width  of  the  diffraction  peak  is  unaffected.  It  is  the  presence 
of  both  positive  and  negative  dislocations  in  the  unit  cell  that  distinguishes  this  case  from 
the  usual  grating  result.  A somewhat  different  physical  insight  derives  from  the  fact  that 
the  scattering  from  adjacent  dislocations  is  phase  shifted  by  7r,  so  forward  scattering  is  zero 
for  qw  = 0,  unlike  the  usual  situation  where  all  elements  of  a structure  scatter  in  phase. 

It  is  necessary  to  avoid  confusion  regarding  the  periodicity  in  g-space  associated  with 
the  wall  of  dislocations,  which  leads  to  a “grating  function”  for  the  dislocations,  and  the 
periodicity  in  the  lattice  which  leads  to  the  Bragg  peaks.  Because  of  the  similarity  in  concept, 
we  will  use  the  term  “dislocation  Bragg”  period  to  refer  to  the  grating  period  of  the  ordered 
(and  partially  ordered)  dislocation  scattering. 

In  the  literature  on  small-angle  scattering  from  particles18,14,17,  Debye  has  given  a simple 
treatment  of  pair  correlation  which  leads  to  analytic  results,  which  is  instructive  to  invoke 
for  the  dislocation  case.  In  particular,  if  it  is  assumed  that  V(r)  has  the  form 

V(r)=  0,  0 < r < Z0, 

V{r)=  1/L2,  l0<r<L , (41) 


then  for  such  a distribution  of  simple  dipoles  or  single-sign  dislocations  in  a single  slip- 
system,  Eqn.  (20)  takes  the  form 


(W)  = blNfit)  + b 2eN(N  - 1 )/2(J)^  /V  - r)  cos  (qwr)  dr 


^JV/2(?)  + ^iV(iV-l)/2(J) 


sin2(qwL/2) 

qlL2/ 4 


— 2(1  — l0/L) 


sin  (qwlo) 


qwL 


(42) 


The  first  term  in  the  brackets  is  the  configuration  factor  for  the  random  wall  (see  Eqn. 
(27)),  and  the  second  term  is  due  to  the  lower  cut-off  at  the  exclusion  distance,  /o-  In 
the  case  of  particles,  the  exclusion  region  exists  because  two  particles  cannot  occupy  the 
same  space.  In  the  case  of  dislocations,  a similar  result  arises  due  to  the  strong  repulsion 
between  two  dislocations  of  like  sign  at  very  small  distances.  For  unlike  sign,  the  dislocations 
will  annihilate  each  other  below  a critical  distance,  even  if  spontaneous  climb  is  necessary. 
This  effect  is  observed  as  a maximum  absolute  density  of  dislocations  in  heavily  deformed 
metals21.  When  2Nl0/L  > 1,  the  configuration  factor  becomes  negative  outside  the  main 
peak,  but  this  can  never  be  the  case,  because  it  would  mean  that  the  mean  distance  between 
dislocations  is  smaller  than  the  exclusion  distance.  For  all  physically-realistic  parameters  in 
Eqn.  (42),  the  effect  of  a small  exclusion  zone  around  a dislocation  is  negligible. 

Next,  consider  the  important  effect  of  partial  ordering  of  defects  in  the  wall.  Rewriting 
Eqn.  (22)  for  a single  slip-system  gives  the  symmetrical  result 


m 

blP(Z)N 
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(43) 


By  symmetry,  we  can  place  two  conditions  on  this  equation 


(w++) 

(w++) 


(w-> 

(W-) 


B= 0 
B——N 


(44) 


Also,  the  joint  pair  distribution  function  for  partially-ordered  dislocations  is  likely  to  re- 
semble that  for  a partially-ordered  two-component  solid.  Thus,  for  dislocations  of  the  same 
sign,  we  consider  the  functional  forms 


V++{r)  = i(l  + Ce-r/Qcos2(7riV+r/L)), 

V—{r)  = -^(l  + Ce-r/acos2{irN-r/L)),  (45) 

where  C is  a parameter  controlling  the  strength  of  the  correlation  and  a is  its  range.  For 
simplicity,  we  do  not  include  an  exclusion  zone  at  the  origin;  such  a zone  would  not  sub- 
stantially affect  the  wall  configuration  factor.  Also,  as  long  as  a « L,  the  error  in  the 
normalization  for  these  equations  will  be  negligible.  The  periods  in  these  functions  were 
chosen  so  that  they  would  be  correct  when  their  contributions  are  significant.  Thus,  when 
B = N,  the  average  spacing  for  the  positive  dislocations  will  be  N/L  = N+/L,  and  when 
B = 0,  the  average  spacing  will  be  N/2L  = N+/L.  Eqn.  (45)  may  become  invalid  for  the 
minority  population  which  will  soon  become  very  highly  disordered,  but  the  contribution  of 
this  term  is  minimized  by  the  pre-factor  (1  — B/N )2  in  (43).  Of  course,  in  these  equations, 
both  C and  a should  be  functions  of  B.  However,  since  we  are  only  interested  in  qualitative 
behavior,  we  will  let  them  be  constants. 

It  remains  to  specify  the  correlation  function  for  unlike  pairs.  The  pre-factor  for  (W+~) 
has  a maximum  at  B = 0,  making  this  the  most  important  point  for  this  term.  Here,  the 
energetics  will  favor  pairing  of  unlike  dislocations  so  that  negative  dislocations  tend  to  sit 
between  positives,  etc.  Thus,  for  simplicity,  we  take 


V + (r)  = -^2  (l  + Ce  T/Q  sin2 (n Nr/ 2L)y  (46) 

In  reality,  the  correlation  strength  and  range  for  unlike  pairs  will  certainly  not  be  the  same 
as  for  like  pairs,  but  such  an  approximation  is  adequate  for  our  purpose  of  studying  the 
qualitative  features  of  partial  ordering.  As  we  move  away  from  B = 0,  the  period  given  in 
Eqn.  (46)  will  remain  correct  because  this  is  the  smallest  period  possible  for  (-1 — ) pairs.  Of 
course,  the  functional  form  may  change  once  a single  sign  dominates,  but  then  the  pre-factor 
for  (W+_)  becomes  very  small,  once  again  minimizing  any  error. 

Fig.  4 shows  the  results  for  this  simplified  version  of  partial  order  near  B = 0.  Here, 
we  plot  only  the  (Wmn)  part  of  the  total  (W)  / b\f2{£)N , i.e.  the  term  proportional  to  N/4 
in  Eqn.  (43).  Figs.  4(a)  and  4(b)  use  B = 0,  N = 20,  and  C = 0.5;  the  other  parameters 
used  are:  (a)  a = 10,  L = 100  and  (b)  a = 20,  L = 200.  As  expected,  the  position  of 
the  dislocation  Bragg  peak  is  set  by  the  dislocation  spacing,  the  width  is  determined  by 
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the  correlation  range,  and  the  strength  depends  on  the  number  of  dislocations  within  this 
range.  Thus,  in  Fig.  4(b),  the  dislocation  Bragg  peak  shifts  position  and  is  narrower,  both 
consistent  with  doubling  the  dislocation  spacing.  Figure  4(c)  uses  the  same  parameters  as 
Fig.  4(a)  except  that  B = N/ 10,  and  shows  the  distinctive  signature  for  the  B ^ 0 case. 
The  narrow  peak  at  the  origin  is  proportional  to  B2  and  it’s  width  is  proportional  to  the 
inverse  of  the  wall  length,  1/L.  This  is  the  same  peak  that  is  found  for  the  ordered  and 
random  cases,  Eqns.  (39)  and  (27).  That  is,  the  peak  at  the  origin  is  independent  of  order, 
and  depends  only  on  the  nonzero  value  of  B. 


IV.  COMPUTER  SIMULATIONS. 

In  this  section,  computer  simulations  are  used  to  confirm  the  correctness  of  the  analysis 
developed  in  the  previous  sections  given  the  various  assumptions  that  were  made.  Cal- 
culations are  performed  on  a square  X-Y  lattice  with  unit  spacing.  This  lattice  does  not 
represent  the  atomic  positions  in  the  sample;  instead,  we  assume  that  a site  of  the  lattice 
represents  a region  of  underlying  crystal.  At  each  site,  we  use  Eqn.  (5)  to  compute  the 
dilatation  of  a distribution  of  dislocations.  Thus,  the  simulation  averages  the  dilatation  over 
a small  region  represented  by  the  lattice  site.  The  resulting  sampling  errors  are  minimized 
by  considering  only  dislocation  distributions  that  are  sparse  relative  to  the  lattice  spacing. 
Since  the  dilatation  is  singular  at  the  core  of  a dislocation,  we  arbitrarily  set  the  dilatation 
of  a cell  containing  a dislocation  center  to  zero.  Eqns.  (1)  and  (2)  are  then  used  to  compute 
the  scattered  intensity  directly  from  the  dilatation  of  the  lattice. 

Three  simulations  are  presented.  In  the  first,  the  form  factor  of  a single  dislocation  is 
examined.  This  is  followed  by  a look  at  the  configuration  factor  of  a random  wall  of  zero 
thickness  with  B = N.  Finally,  the  case  of  a random  wall  of  zero  thickness  with  B = 0 is 
investigated.  The  random  wall  results  are  obtained  by  using  a random  number  generator  to 
generate  a large  number  of  random  dislocation  configurations  with  the  required  parameters. 
The  configuration  factors  calculated  from  these  configurations  are  then  averaged  to  obtain 
the  final  results.  In  all  of  the  simulations,  the  Burgers  vectors  of  the  dislocations  have 
equal  magnitude  and  lie  in  the  Y-direction.  Dislocation  walls  are  aligned  along  the  Y- 
axis,  perpendicular  to  the  slip  plane.  In  this  configuration,  the  qw  direction  is  favorable  for 
scattering,  but  no  scattering  will  take  place  in  the  thickness  direction  because  of  the  single 
dislocation  form  factor. 

The  first  simulation  is  meant  to  confirm  the  scattering  from  a single  dislocation  as  given 
by  Eqn.  (7).  The  dislocation  is  located  at  the  center  of  a square  lattice  with  an  edge  length 
of  400;  the  results  are  shown  in  Fig.  5.  The  units  for  qw  are  27t /d,  where  d is  a length 
measured  in  units  of  lattice  spacings.  First,  note  that  the  scattered  intensity  goes  to  zero 
as  qw  — > 0,  and  the  expected  oscillatory  Bessel  function  character  is  confirmed.  Also,  the 
maximum  of  the  scattering  occurs  at  qw  = 2n/R  as  predicted  by  Eqn.  (7)  with  the  upper 
cut-off  distance,  R,  set  to  the  distance  from  the  dislocation  to  the  sample  edge.  Finally, 
the  log  plot  shown  in  Fig.  5(b)  shows  that  the  asymptotic  form  of  the  scattering  is  l/q^, 
for  large  qw:  as  predicted.  Not  shown  is  the  scattering  for  qt , which  is  strictly  zero,  as 
predicted  by  Eqn.  (7).  Although  the  general  behavior  of  the  computer  simulation  matches 
the  behavior  predicted  by  Eqn.  (7),  we  note  that  this  agreement  is  not  quantitative.  The 
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reason  for  the  difference  is  that  Eqn.  (7)  was  derived  with  the  assumption  that  finite-size 
effects  are  negligible.  This  is  clearly  not  true  for  the  computer  simulations. 

It  is  important  to  note  that  the  oscillations  predicted  by  Eqn.  (7)  and  shown  in  Fig.  5 
would  not  be  observed  in  an  actual  experiment.  This  is  because  the  coherence  length,  R,  is 
not  a sharp  cut-off  that  is  the  same  for  all  dislocations.  Thus,  averaging  over  the  dislocations 
in  a sample  would  produce  an  “averaged”  structure  factor  with  a simple  q~ 2 behavior  and 
no  oscillations. 

It  is  useful  to  examine  what  happens  to  the  scattering  at  qw  = 0 when  the  dislocation 
is  displaced  from  the  origin  in  the  Y-direction.  This  is  related  to  the  qw  = 0 scattering  for 
a dislocation  dipole.  When  a single  dislocation  is  displaced  normal  to  its  Burgers  vector 
in  the  T-direction,  since  the  computation  cell  is  finite,  the  net  volume  change  in  the  cell  is 
non-zero.  This  is  because  the  number  of  atoms  above  the  slip  plane,  where  the  dilatation  is 
negative,  is  different  from  the  number  below  the  slip  plane,  where  the  dilatation  is  positive. 
Thus,  the  scattering  at  qw  = 0 increases  as  the  dislocation  is  displaced  in  the  positive  In- 
direction. Similarly,  for  a dislocation  dipole  of  positive  and  negative  Burgers  vector,  with 
the  dipole  oriented  along  the  Y -axis,  there  is  again  a net  volume  change  for  the  computation 
cell,  overall.  Thus,  a wall  of  alternating  dislocations  along  the  Y-axis  will  have  a net  volume 
change  equivalent  to  a single  dipole  of  effective  displacement  equal  to  the  length  of  the  wall 
and  effective  Burgers  vector  equal  to  6/2  (or  equivalently,  a single  whole  dipole  with  a length 
equal  to  half  the  wall  length). 

The  plots  in  Fig.  6 show  simulations  of  the  wall  function,  Eqn.  (9),  for  randomized  walls 
of  zero  thickness  with  100  dislocations,  which  are  averaged  over  1000  independent  iterations. 
In  Fig.  6(a),  all  of  the  dislocations  have  the  same  sign  and  in  Fig.  6(b),  the  dislocations 
have  zero  net  Burgers  vector.  These  simulations  should  be  compared  with  Eqns.  (27)  and 
(26).  In  both  cases,  we  have  set  \be\  = /(£)  = 1.  Fig.  6(a)  for  the  same  sign  wall  shows  that 
the  scattering  has  a strong  peak  at  the  origin  proportional  to  the  square  of  the  number  of 
dislocations.  (Of  course,  the  total  scattered  intensity  is  a product  of  the  single  dislocation 
function,  which  goes  to  zero  at  qw  = 0,  and  otherwise  multiplies  the  wall  function  by  q~ 2 
when  qw  is  large.)  At  large  q,  the  function  approaches  it’s  asymptotic  value  of  N.  The 
B = 0 wall  function,  depicted  in  Fig.  7(b),  has  a “hole”  at  the  origin,  and  then  rises  to 
the  same  asymptotic  value  as  before.  In  both  Figs.  6(a)  and  6(b)  the  peak  (and  hole)  have 
widths  given  by  5qw  = 2tt/L,  as  predicted.  We  note  that  the  wall  function  in  both  cases 
approaches  an  asymptotic  value  proportional  to  the  total  number  of  dislocations  in  the  wall. 
This  asymptotic  value  is  the  term  we  have  called  the  “self  scattering”  term  in  Eqn.  (12) 
and  its  later  incarnations. 


V.  CONCLUSIONS. 

The  theory  developed  in  this  paper  for  small-angle  scattering  by  dislocations  in  deformed 
single-crystals  has  the  traditional  feature  that  the  scattering  is  a product  of  several  inde- 
pendent factors.  Aside  from  the  atomic  structure  factor  and  the  sample  size  form  factor, 
neither  of  which  is  considered  in  this  paper,  the  intensity  is  composed  of  a structure  factor 
for  the  “single  defect” , which  may  be  either  a single  dislocation  or  a dislocation  dipole  pair, 
and  a configuration  factor  for  the  configuration  of  single  defects. 
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Experiments  capable  of  measuring  the  relevant  parameters  will  be  difficult  and  very  few 
existing  facilities  currently  have  the  necessary  equipment.  Experimental  difficulties  include 
such  requirements  as  avoidance  of  accidental  Bragg  peaks,  adequate  rotational  control  of  the 
sample/beam  geometry,  high-quality  single-crystal  dog-bone  shaped  metal  samples  that  are 
thin  enough  to  limit  multiple  scattering,  a tensile  stage  capable  of  deforming  these  samples  in 
situ  while  simultaneously  measuring  the  engineering  stress  and  strain,  high  beam  brightness, 
and  ultra-low-angle  capabilities.  The  smallest  q that  can  be  reached  using  existing  suitable 
USAXS  facilities  corresponds  to  a spatial  distance  of  « 1/im.  Thus,  large  structures  such  as 
entire  dislocation  walls  and  cells  are  not  presently  visible.  However,  the  local  organization  of 
dislocations  within  a wall  should  certainly  be  accessible.  Thus,  the  theory  should  be  applied 
to  the  structure  of  walls,  and  to  the  scattering  of  individual  dislocations,  but  not  to  the 
scattering  by  complete  cells.  A new  USAXS  facility  under  construction  at  the  Advanced 
Photon  Source  will  extend  this  usable  g-range  by  about  a factor  of  5,  which  should  allow 
the  observation  of  complete  cells  in  many  cases. 

The  structure  factor  for  the  single  dislocation  can  be  written  as  two  separate  factors 
corresponding  to  the  uqz  scattering”  for  the  scattering  due  to  the  length  coordinate  of  the 
dislocation,  and  the  ^-scattering  normal  to  the  length  direction.  The  ^-structure  factor  has 
a strong  and  narrow  maximum  when  the  scattering  vector  is  normal  to  the  length  direction, 
with  a simple  “diffraction  function”  angular  dependence.  The  scattered  intensity  reaches  a 
maximum  when  the  scattering  vector  is  also  perpendicular  to  the  edge  components  of  the 
Burgers  vectors.  For  the  single  dislocation,  this  single  defect  angular  variation  is  rather 
“soft”  with  the  form  sin2£,  while  the  dislocation  dipole  is  more  highly  directional  with  a 
“spot  light”  in  the  direction  normal  to  the  slip  plane  of  form  sin4  £.  Thus,  a dislocation  is 
best  situated  to  scatter  when  the  scattering  vector  is  perpendicular  to  its  slip  plane. 

The  single  dislocation  structure  factor  is  multiplied  by  the  factor  for  the  configuration 
of  the  dislocations,  which  we  have  assumed  is  a diffuse  wall.  The  wall  function,  (W),  is 
composed  of  two  additive  terms.  The  first  term  is  the  self-scattering  term  which  goes  to 
zero  as  qwL  0.  In  both  of  the  specific  cases  examined  (ie,  the  random  case  and  the  fully 
ordered  case),  this  term  goes  to  zero  quadratically.  For  qwL  » 1,  this  first  term  in  the  wall 
function  is  independent  of  qw,  and  is  simply  proportional  to  the  total  number  of  scattering 
dislocations  (or  dislocation  dipoles).  Because  it  is  a constant  for  large  qw,  this  term  in 
the  wall  function  will  be  dominant  when  other  contributions  that  decrease  with  qw  have 
decayed.  In  the  large  qw  domain,  then,  the  scattering  is  characterized  simply  by  the  single 
defect  structure  factor,  without  the  complication  of  the  defect  configuration,  and  the  strength 
of  the  scattering  is  linear  in  the  total  number  of  scatterers.  An  important  consequence  of 
this  fact  is  that  if  single  dislocations  and  dipoles  are  both  present  in  significant  numbers, 
the  problem  becomes  solvable  in  this  qw  range.  Thus,  the  scattered  intensity  in  this  range 
will  be  the  sum  of  q~ 2 for  single  dislocations  and  q~3  for  dislocation  dipoles.  Experiments 
in  this  range  could  therefore  settle  the  long-standing  debate  on  the  existence  of  dipoles  in 
deformation  structures.2,7,8 

The  second,  g-dependent,  term  comprising  the  wall  function  is  again  composed  of  two 
separate  factors  that  are  multiplied  together.  The  first  factor  describes  the  dimension  along 
the  wall  (i.e.  the  qw  direction),  and  perpendicular  to  the  dislocation  lines.  (In  the  case  of  fee 
crystals,  it  is  reasonable  to  assume  that  the  walls  are  composed  of  dislocations  parallel  to 
one  another,  but  not  necessarily  on  the  same  slip  planes.)  This  dimension  was  normal  to  the 
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slip  plane  in  the  computer  simulations.  This  factor  depends  both  on  the  net  Burgers  vector 
of  the  wall  and  on  the  presence  of  ordering  in  the  qw  direction.  For  a completely  random 
distribution  with  B = N,  this  function  has  a peak  at  the  origin  proportional  to  B 2,  which 
decreases  quadratically  to  zero  as  a function  of  qw  in  the  range  qwL  1.  For  a completely 
random  distribution  with  B = 0,  the  term  is  zero. 

When  the  dislocations  are  partially  ordered,  interference  peaks  appear  associated  with 
the  inter-dislocation  distance  and  the  correlation  length.  For  B ^ 0,  the  same  peak  at  the 
origin  appears  which  was  discussed  for  the  random  case.  For  all  B,  a dislocation  Bragg 
peak  exists  at  a position  characterized  by  the  defect  structure  repeat  distance.  For  complete 
ordering,  the  peak  is  sharp  with  a width  given  by  the  number  of  dislocations  in  the  wall, 
but  if  the  order  is  only  short  range,  with  correlations  extending  only  a few  inter-dislocation 
distances,  then  the  peak  is  much  broader. 

The  second  factor  in  the  ^-dependent  term  describes  the  ^-dependence  normal  to  the 
wall  in  the  thickness  direction.  The  thickness  configuration  factor  is  a delta  function  when 
the  thickness  goes  to  zero,  and  otherwise  depends  on  the  distribution  of  the  dislocations  in 
the  thickness  direction.  If  the  dislocations  are  distributed  randomly  within  a thick  wall,  and 
the  wall  has  a sharp  interface,  then  the  qt  scattering  will  have  a conventional  q t~2  behavior. 
If  the  interface  is  diffuse  rather  than  sharp,  the  power  will  deviate  from  — 2 to  — (2  + £).  5 
contains  information  on  the  width  of  the  diffuse  boundary. 

In  summary,  any  experimental  investigation  of  the  dislocation  structures  must  be  guided 
by  geometrical  considerations  of  the  incident  X-ray  or  neutron  direction,  and  its  relation  to 
the  relevant  slip  systems  of  the  material.  It  is  found  that  only  those  dislocations  lying  normal 
to  the  scattering  vector  will  contribute  significantly  to  the  scattering  due  to  the  strong  peak 
in  the  configuration  factor  for  the  single  dislocation.  Further,  the  scattering  will  be  peaked 
(very  strongly  in  the  case  of  dislocation  dipoles)  for  q also  perpendicular  to  the  Burgers 
vector.  Thus,  the  most  favorable  scattering  geometry  will  be  one  with  the  scattering  vector 
perpendicular  to  the  slip  plane  of  the  dominant  dislocations.  Above  a particular  |q|,  the 
scattering  will  be  sensitive  only  to  the  total  number  of  edge  dislocations,  while  below  this 
value,  the  scattering  will  be  sensitive  to  the  ordering  within  the  walls.  For  carpet  structures, 
the  scattering  will  be  determined  by  the  distribution  in  the  thickness  direction,  while  for  walls 
normal  to  the  dominant  slip  plane,  the  scattering  will  be  determined  by  the  ordering  within 
the  plane  of  the  walls  and  the  net  Burgers  content.  Structural  features  will  be  accessible 
only  to  the  extent  that  the  lowest  measurable  |q|  corresponds  to  the  length  scale  of  the 
relevant  structure.  Thus,  dislocation  structure  within  a wall  should  be  readily  accessible  by 
current  USAXS,  while  the  cellular  structure  formed  by  the  walls  will  be  accessible  only  by 
the  next  generation  facilities  and  for  specimens  at  relatively  high  strain  levels.  In  theory, 
it  should  be  possible  to  obtain  considerable  information  about  dislocation  structures  in  a 
deformed  metal  by  a considered  variation  of  the  geometry  between  the  incident  beam  and 
the  active  slip  systems  of  the  material.  It  remains  to  be  seen  just  how  much  information 
can  be  extracted  from  such  experiments  in  practice.  This  question  is  being  pursued  further. 
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FIGURES 


FIG.  1.  Coordinate  system  for  dislocation.  The  dislocation  lies  along  the  Z-axis  perpendicular 
to  the  paper  at  the  origin,  and  the  edge  component  of  the  Burgers  vector  lies  along  the  X-axis. 
The  projection  of  the  scattering  vector,  q,  in  the  X-Y  plane  has  cylindrical  coordinates  (<?p,£).  A 
point  in  the  X-Y  plane  has  the  cylindrical  coordinates,  (r,  6).  The  wall  extends  in  the  w direction 
with  t describing  the  perpendicular  direction. 


FIG.  2.  A dislocation  dipole  at  the  origin  is  oriented  along  the  Y-axis,  with  separation  distance 
2d.  The  edge  component  of  the  Burgers  vectors  of  the  dipole  is  along  the  X-axis.  The  coordinate 
definitions  are  the  same  as  those  in  Fig.  1. 


FIG.  3.  Coordinate  system  for  integration  over  wall  variables.  In  (a),  the  original  integration 
region  in  the  variables  (R,Rr)  is  a square  of  side  L.  In  the  new  variables,  (r,  r'),  the  integration 
region  is  the  rotated  rhombus,  (b). 


FIG.  4.  Partially  ordered  wall,  (a)  Plot  of  term  proportional  to  iV/4  in  Eqn.  (42)  with  L = 100, 
C = 0.5,  q = 10,  N = 20,  and  B = 0.  The  plot  shows  the  dislocation  Bragg  peak  with  the  position 
set  by  the  dislocation  spacing,  the  width  set  by  the  correlation  distance,  and  the  strength  by  the 
number  of  correlated  dislocations,  (b)  Same  as  above,  with  L — 200  and  a = 20.  The  peak  has 
shifted  position  to  half  its  previous  -value.  Its  peak  height  is  the  same,  because  the  same  number 
of  correlated  dislocations  are  present  in  each  case  and  so  is  the  total  number  of  dislocations  in  the 
wall.  The  peak  width  is  reduced  by  half,  because  the  correlation  distance  is  increased,  (c)  Same 
as  (a),  with  B/N  = 0.1.  Now  a small  narrow  peak  appears  at  the  origin.  As  B is  increased,  this 
peak  grows  proportionally  to  B 2. 


FIG.  5.  Scattering  for  a single  dislocation  at  the  origin  of  a square  lattice  with  401  sites  on  a 
side,  a)  is  a cut  along  the  g^-axis,  showing  the  intensity  of  scattering  as  a function  of  q.  b)  is  a 
log  plot  of  the  same  data  showing  that  the  qw  dependence  for  large  qw  is  q~2. 


FIG.  6.  Wall  configuration  factors  for  random  walls  containing  100  dislocations  with  (a)  B = N 
and  (b)  B = 0.  The  figures  show  (VV)  plotted  against  qw.  The  wall  length  is  500  lattice  spacings. 
Dislocation  positions  are  randomized  and  averaged  over  1000  simulations.  In  both  cases,  the 
asymptotic  value  for  large  qw  is  the  total  number  of  dislocations. 
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